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ABSTRACT 

Objective: To examine the effects of designation and 
volume of neonatal care at the hospital of birth on 
mortality and morbidity outcomes in very preterm 
infants in a managed clinical network setting. 
Design: A retrospective, population-based analysis of 
operational clinical data using adjusted logistic 
regression and instrumental variables (IV) analyses. 
Setting: 165 National Health Service neonatal units in 
England contributing data to the National Neonatal 
Research Database at the Neonatal Data Analysis Unit 
and participating in the Neonatal Economic, Staffing 
and Clinical Outcomes Project. 
Participants: 20 554 infants born at <33 weeks 
completed gestation (17 995 born at 27-32 weeks; 
2559 born at <27 weeks), admitted to neonatal care 
and either discharged or died, over the period 1 
January 2009-31 December 201 1 . 
Intervention: Tertiary designation or high-volume 
neonatal care at the hospital of birth. 
Outcomes: Neonatal mortality, any in-hospital 
mortality, surgery for necrotising enterocolitis, surgery 
for retinopathy of prematurity, bronchopulmonary 
dysplasia and postmenstrual age at discharge. 
Results: Infants born at <33 weeks gestation and 
admitted to a high-volume neonatal unit at the hospital 
of birth were at reduced odds of neonatal mortality (IV 
regression odds ratio (OR) 0.70, 95% CI 0.53 to 0.92) 
and any in-hospital mortality (IV regression OR 0.68, 
95% CI 0.54 to 0.85). The effect of volume on any in- 
hospital mortality was most acute among infants born 
at <27 weeks gestation (IV regression OR 0.51, 95% CI 
0.33 to 0.79). A negative association between tertiary- 
level unit designation and mortality was also observed 
with adjusted logistic regression for infants born at 
<27 weeks gestation. 

Conclusions: High-volume neonatal care provided at 
the hospital of birth may protect against in-hospital 
mortality in very preterm infants. Future developments of 
neonatal services should promote delivery of very preterm 
infants at hospitals with high-volume neonatal units. 



Strengths and limitations of this study 



A national dataset consisting of the electronic 
patient records of a large majority of admissions 
to neonatal specialist care in England. 
The analysis takes into account observed and 
unobserved confounding. 
A weakness is that the analysis is unable to dis- 
entangle the effects of the neonatal unit at the 
place of birth from subsequent transfers to other 
neonatal units. 



INTRODUCTION 

Intense debate has revolved around the 
optimal organisation of neonatal critical care 
services. Numerous studies have suggested 
that the intensity and volume of neonatal care 
at the hospital of birth is negatively correlated 
with adverse clinical outcomes, including 
mortality. 1-11 This has contributed to calls for 
centralisation of neonatal services and the 
closure of smaller neonatal units. 2 11 12 

Following a review by the Department of 
Health in 2003, perinatal centres in England 
were reorganised into managed clinical net- 
works (MCN). 13 MCNs provide some of the 
benefits of centralisation, but also strive to 
maintain equity and ease of access to services 
by keeping lower care level and lower volume 
neonatal units open, with provision for trans- 
fer to higher care level or higher volume 
units, if required. 13 Particular emphasis is 
placed on the importance of transferring 
women at risk of extremely preterm labour 
to tertiary centres before delivery. 
Consequently, most networks aim to transfer 
women at high risk of delivery at <27 weeks 
gestation. We have previously shown that, 
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since the formation of MCNs, the proportion of low- 
gestational age infants born in hospitals with higher des- 
ignation neonatal units and their transfer rate between 
hospitals has increased significantly; however, it remains 
unclear what effect this has had on clinical outcomes. 14 

Studies that have examined the effects of neonatal 
unit designation or volume of neonatal care provided at 
the hospital of birth have shown that low designation 
level or volume is associated with increased rates of mor- 
tality, 1-10 decreased infection rate, 7 increased severe 
periventricular haemorrhage 11 and increased broncho- 
pulmonary dysplasia. 7 However, these studies were 
almost exclusively conducted in the USA where there is 
greater variability in neonatal unit volume — the highest 
volume units in the USA are typically much larger than 
equivalent units in England — and there are no formal 
arrangements for MCNs. Results from similar studies 
using data from the UK are limited and based on data 
from 1998 to 1999, prior to the formation of MCNs. 15 16 
We are not aware of any studies that have examined 
infant outcomes for neonatal specialist services in MCNs 
in relation to unit designation or volume. In addition, 
organisation of neonatal care differs between countries 
potentially affecting the generalisability of results from 
these systems; for example, in Germany neonatal ser- 
vices are markedly deregionalised whereas in Finland 
and Portugal there is a high degree of regionalisation. 17 

Our aim in this study was to examine the effects of 
designation and volume of neonatal care provided at the 
hospital of birth on mortality and morbidity outcomes. 
We assess whether organisational factors remain determi- 
nants of clinical outcomes despite the goals of neonatal 
reorganisation that sought to ensure that vulnerable 
infants are not disadvantaged by their place of birth. 



METHODS 

Data source and study population 

For the purpose of this empirical investigation, we 
extracted data from the National Neonatal Research 
Database (NNRD) for neonatal units participating in the 
Neonatal Economic, Staffing and Clinical Outcomes 
Project (NESCOP). The NNRD is held by the Neonatal 
Data Analysis Unit (NDAU), Imperial College, London, 
and was created from patient-level electronic records of 
all infants admitted to 168 of 173 neonatal units in 
England. NESCOP included 165 centres providing peri- 
natal care. On behalf of NESCOP, the Medical Research 
Council (MRC) EPICure studies carried out the Unit 
Profile Survey (UPS) during 2011, comprising a survey 
of English hospitals that provided onsite obstetric and 
neonatal services. We extracted records from the NNRD 
of all infants born in participating centres at <32 +6 
weeks +days gestation, admitted over the period 1 January 
2009-31 December 2011, and who were discharged or 
died over the same period. We excluded infants who 
only received transitional care (n=5), which was defined 



according to English Department of Health's Healthcare 
Resource Group (HRG4) code 'XA04Z'. 18 Gestational 
age was determined by ultrasound scan. 

Outcomes 

We derived the following outcomes from the extracted 
data for use in the analyses: 28-day (neonatal) mortality, 
any in-hospital mortality, surgery for necrotising entero- 
colitis, treatment for retinopathy of prematurity (ROP) 
and bronchopulmonary dysplasia (BPD). We defined 
BPD as the requirement of supplementary oxygen for at 
least 28 days and at 36 weeks postmenstrual age 
(PMA). 19 We also examined PMA at discharge as a 
marker of length of stay; this was defined as the gesta- 
tional age at birth plus the length of stay at final dis- 
charge from any neonatal unit or death. We defined the 
outcome to be one if the PMA at discharge was greater 
than 40 weeks and zero otherwise. 

Covariates 

To determine appropriate covariates, we reviewed previ- 
ous prediction models for very preterm infants 20 and 
selected variables that (1) were significant predictors of 
adverse sequelae, (2) were available in our dataset and 
of high quality and (3) not confounded by the provision 
of neonatal care. The variables we included were: gesta- 
tional age at birth, gestational age squared, birthweight 
z-score (birth weight standardised by gestational age 
week) and the following indicators: whether the mother 
received a full or partial course of antenatal steroids, 
sex, infant year of birth and whether or not the mother 
came from an area within the lowest decile of the Index 
of Multiple Deprivation 2007 score. 21 

Statistical methods 

We conducted two separate sets of analyses based on 
whether or not infants were admitted to a neonatal unit 
at the hospital of birth designated as: (1) a tertiary 
centre 22 or (2) high volume. For the latter, we defined 
volume according to the annual number of care days at 
any level of care provided to very preterm infants (<32 +6 
weeks gestation). A 'high-volume' unit was defined as 
one whose volume was in the top quartile of all neonatal 
units in the sample. 'High volume' was determined by 
quartile rather than an absolute care day threshold to 
facilitate comparison with other measures of volume in 
the sensitivity analyses. A previous study that examined 
organisational characteristics of neonatal units also cate- 
gorised volume using quartiles. 17 Dichotomising by 
upper quartile divided the infants between high-volume 
and low-volume units in approximately the same propor- 
tion as between tertiary-level and non-tertiary-level units. 
To aid comparison with other studies, in particular from 
the USA, and as a robustness check, 'high volume' was 
also defined as 100 very low birthweight (VLBW; 
<1500 g) admissions of infants born in the same hospital 
per annum. 
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We first conducted an unadjusted comparison of clin- 
ical characteristics and outcomes of infants by unit 
characteristics. Second, we estimated an adjusted model 
and finally, we conducted an adjusted comparison using 
an instrumental variables methodology to account for 
unobserved confounding. In the absence of a rando- 
mised control trial, instrumental variables methodology 
acts as an ex post randomisation and enables us to esti- 
mate the 'causal effects' of designation and volume of 
neonatal care provided at the hospital of birth. The 
methodology involves the use of a variable called an 
'instrument' which, in this context, needs to fulfil two 
criteria: (1) it should be strongly correlated with the 
characteristics of the neonatal unit at the hospital of 
birth and (2) it should be uncorrelated with the out- 
comes of interest conditional on observed covariates and 
therefore uncorrelated with unobserved confounders. 

For the instruments, we used indicators for the desig- 
nated level of care of the nearest neonatal unit to the 
mother's residence, an indicator for whether it had sur- 
gical facilities, an indicator for whether it was high 
volume, the distance to the nearest neonatal unit and 
the interactions of either the level of care indicators or 
high-volume indicator with distance, giving nine instru- 
ments in total. Straight line distance was calculated from 
the population-weighted centre of the mother's lower 
super output area to each hospital. 23 

These instrumental variables fulfil condition (1) if 
infants are more likely to be born in the hospital closest 
to the mother's residence. They will also fulfil condition 
(2) if the location of the mother's residence is uncorre- 
lated with an infant's unobserved clinical risk. We tested 
for a difference in observed characteristics by level and 
volume of the nearest neonatal unit. However, tertiary- 
level and high-volume units are more likely to be in 
urban areas that are socioeconomically deprived so we 
may expect to see more preterm and low birthweight 
infants being born in these areas. 24 We therefore also 
controlled for local deprivation when testing for a differ- 
ence in means by nearest neonatal unit characteristics 
by estimating a linear regression of the observed variable 
of interest on the nearest neonatal unit characteristic 
and deprivation indicator, and using an F-test to test the 
coefficient on the nearest neonatal unit characteristic 
variable. 

As the outcomes are all binary logistic regression was 
used. In order to employ instrumental variables estima- 
tion in this framework, two-stage residual inclusion 
(2SRI) was used. 25 The 2SRI method is explained in 
online supplementary appendix A. The SEs were 
adjusted for clustering within units. 

Our baseline analyses examined infants born at <32 +6 
weeks gestation. We then conducted analyses on subsets 
of infants born at <26 +6 weeks gestation or at 27 +0 -32 +6 
weeks gestation; <26 +6 weeks gestation is the cut-off used 
by perinatal networks for prioritising interunit transfers. 
'Statistical significance', where discussed, refers to a 5% 
significance level in all cases. 



Missing data and sensitivity analyses 

Infants with missing outcomes data were excluded from 
the analyses, while those with missing covariate data 
were assigned a zero in the case of binary indicators. 
There were no infants with missing continuous covari- 
ates. We excluded all infants with any missing data as a 
further sensitivity analysis. 

Separate sensitivity analyses, using our preferred 
method of instrumental variables logistic regression, also 
explored the effects of: (1) including unit random 
effects in the statistical models; (2) removing infants 
who died from analyses of the morbidity and PMA at dis- 
charge outcomes and defining a new outcome of any 
in-hospital mortality and/ or BPD to account for possible 
bias caused by infants dying prior to experiencing the 
morbidity outcome; (3) redefining high volume as the 
top 25% of units in terms of intensive care days provided 
to <32 +6 gestational week infants; (4) redefining high 
volume as the top 25% of units in terms of number of 
<32 +6 gestational week infants cared for and (5) redefin- 
ing high volume as at least 100 VLBW infants born in 
and admitted to the neonatal unit in the hospital per 
annum. 

All analyses were carried out with R V.2.14.2 and Stata 
V.ll. 



RESULTS 

In total, data for 20 554 infants born at <32 +6 weeks ges- 
tation over the study period and admitted to a neonatal 
unit at the hospital of birth were extracted from the 
NNRD, 2559 of whom were born at <26 +6 weeks gesta- 
tion. Table 1 provides descriptive statistics of the samples 
analysed. 

In the sample, 9466 (46.1%) infants were born in hos- 
pitals with a tertiary-level neonatal unit and 9541 
(46.4%) were born in hospitals with a high-volume neo- 
natal unit. The cut-off for high volume was approxi- 
mately 3480 annual care days for infants born at <32 +6 
weeks gestation in each hospital. The total sample of 
20 554 infants were born in 165 different hospitals, 44 
(26.7%) of which had level 3 neonatal units, 81 (49%) 
level 2 neonatal units and 39 (23.6%) level 1 neonatal 
units. There were 39 (23.6%) neonatal units classified as 
high volume, 30 (78%) of which were designated level 3 
units; consequently, 14 of the 44 (31.8%) level 3 desig- 
nated units were not classified as high volume. Among 
the 20 554 infants, 1892 (9.2%) were born in hospitals 
with neonatal units that were classified as high volume 
but not tertiary level and 1817 (8.8%) were born in hos- 
pitals with neonatal units classified as tertiary level but 
not high volume. 

'Standard' adjusted results 

Table 2 presents the estimated adjusted ORs associated 
with admission to either tertiary or high-volume neo- 
natal care at the hospital of birth. 
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Table 1 Descriptive statistics for preterm infants born <32 +6 weeks gestation by neonatal unit characteristic at the hospital 
of birth 



Designation of unit 



Volume of unit* 



Tertiary-level 
unit 



Non-tertiary- 
level unit 



p Valuet 



High- 
volume unit 



Non-high- 
volume unit 



p Valuef 



n(%) 9466 (46.1) 

Gestation (weeks), mean (SD) 29.2 (2.5) 

Birth weight (g), mean (SD) 1313.9 (438.7) 

Received full or partial course of 6394 (67.6) 
antenatal steroids 

Deprivation score bottom 10% 2020 (21 .4) 

Male 5048 (53.3) 

Neonatal mortality 423 (4.5) 

Any in-hospital mortality 569 (6.0) 

BPDJ 3695 (39.0) 

Treatment for ROP 226 (2.4) 

Surgery for NEC 167(1.8) 

PMA§ at discharge >40 +0 weeks 1 292 (1 3.7) 



1 1 088 (54.0) 

30.0(2.1) <0.001 

1451.9 (404.5) <0.001 

7262 (65.5) 0.002 



9541 (46.4) 11 013 (53.6) 

29.3 (2.5) 29.9 (2.2) <0.001 

1 326.6 (436.7) 1 441 .8 (409.4) <0.001 

6330 (66.4) 7326 (66.5) 0.790 



1342 (12.1) 
5397 (53.4) 
366 (3.3) 
425 (3.8) 
2856 (25.8) 
107 (1.0) 
123 (1.1) 
848 (7.7) 



<0.001 
0.756 
<0.001 
<0.001 
<0.001 
<0.001 
<0.001 
<0.001 



1730 (18.1) 

5093 (53.4) 
394 (4.1) 
527 (5.5) 

3548 (37.2) 
195 (2.0) 
163 (1.7) 

1237 (13.0) 



1632 (14.8) 
5892 (53.5) 
395 (3.6) 
467 (4.2) 
3003 (27.3) 
138 (1.3) 
127 (1.2) 
903 (8.2) 



<0.001 
0.863 
0.043 
<0.001 
<0.001 
<0.001 
0.001 
<0.001 



All values are n (%) unless otherwise stated. 

*High volume was defined as being in the top quartile of units by number of care days provided to infants born at <32 +6 weeks gestation. 
fContinuous variables were tested by t test, categorical variables by % 2 test. 

jiBPD defined as requirement of supplementary oxygen for at least 28 days postbirth and at 36 weeks postmenstrual age. 
§PMA at discharge equal to gestational age at birth plus length of stay in weeks. 

BPD, bronchopulmonary dysplasia; NEC, necrotising enterocolitis; PMA, postmenstrual age; ROP, retinopathy of prematurity. 



The standard logistic regressions did not reveal a statis- 
tically significant difference in the OR of mortality for 
very preterm infants admitted to tertiary-level care at the 
hospital of birth compared with their counterparts 
admitted to non-tertiary-level care. However, when con- 
sidering only infants born at <26 +6 weeks gestation, we 
found a reduction in the OR of neonatal mortality (OR 
0.65, 95% CI 0.46 to 0.91, p=0.012), but not any 
in-hospital mortality. 

For infants admitted to a high-volume neonatal unit at 
the hospital of birth, a reduced OR of neonatal mortality 
was observed for those born at <32 +6 weeks gestation 
(OR 0.73, 95% CI 0.56 to 0.95, p=0.018) and at <26 +6 
weeks gestation (OR 0.62, 95% CI 0.44 to 0.87, 
p=0.006), but this was not replicated for infants born at 
27 +0 to 32 +6 weeks gestation. Those infants born at <26 +6 
weeks gestation were also at reduced OR of any 
in-hospital mortality (0.71, 95% CI 0.52 to 0.97, 
p=0.033) and increased OR of BPD (OR 1.59, 95% CI 
1.18 to 2.14, p=0.002) compared with their counterparts 
admitted to a non-high-volume neonatal unit at the hos- 
pital of birth. There were no other statistically significant 
differences observed for the morbidity outcomes. 

Instrument validity 

The instruments were strongly correlated with the 
characteristics of the unit at the hospital of birth; 88.4% 
of infants whose nearest neonatal unit was designated 
level 3 were born in a hospital with a level 3 unit com- 
pared with only 22.5% of infants whose nearest neonatal 
unit was not designated level 3. Table 3 shows descriptive 
statistics for the 20 554 very preterm infants by the 



designation and volume of the neonatal unit nearest to 
the mother's place of residence. After correcting for 
deprivation, there were no statistically significant differ- 
ences in the observed covariates. 

Instrumental variables logistic regression 

Table 4 shows the estimated ORs using the instrumental 
variables logistic regressions. We found no significant dif- 
ferences in neonatal mortality between infants admitted 
to either tertiary or non-tertiary neonatal care at the hos- 
pital of birth. We did find an increased OR of treatment 
for ROP for very preterm infants born at 27 +0 -32 +6 
weeks gestation born in a hospital with a tertiary-level 
unit (OR 2.17, 95% CI 1.06 to 4.47, p=0.035). 

In contrast to the effect of tertiary-level care, admis- 
sion to a high-volume neonatal unit at the hospital of 
birth significantly reduced the OR of neonatal mortality 
(OR 0.70, 95% CI 0.53 to 0.92, p=0.011) and any 
in-hospital mortality (OR 0.68, 95% CI 0.54 to 0.85, 
p=0.001) in very preterm infants. These effects were 
most acute among infants born at <26 +6 weeks gestation. 
In terms of morbidity, the only significant effect was 
found for BPD (OR 1.78, 95% CI 1.12 to 2.81, p=0.014) 
for infants born at <26 +6 weeks gestation and admitted 
to high-volume neonatal care at the hospital of birth. 

Sensitivity analyses 

The results from the sensitivity analyses are presented in 
online supplementary appendix B. There were 1172 
(5.7%) infants with missing data for antenatal steroids; 
there were no missing values for the other covariates. 
The results remained qualitatively similar when all 
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Table 2 Adjusted ORs for outcomes associated with admission to either tertiary or high-volume neonatal care at the hospital of birth using a 'standard' logistic regression 
model 



Tertiary neonatal unit High-volume neonatal unitt 



Outcome 


(1) 

<32 +6 weeks 


(2) 

<26 +6 weeks 


(3) 

27+0 _ 32+6 weeks 


(4) 

<32 +6 weeks 


(5) 

<26 +6 weeks 


(6) 

27+0 _ 32+6 weeks 


Neonatal mortality 
Any in-hospital mortality 
BPD 

Treatment for ROP 

Surgery for NEC 

PMA at discharge >40 weeks 


0.77 (0.59 to 1 .00) 
0.91 (0.72 to 1.15) 
1 .23** (1 .07 to 1 .40) 
1 .26 (0.91 to 1 .75) 
1 .05 (0.76 to 1 .44) 
1.17 (0.97 to 1.41) 


0.65* (0.46 to 0.91) 
0.78 (0.57 to 1 .06) 

1.50** (1.11 to 2.01) 
1 .09 (0.76 to 1 .57) 
0.89 (0.58 to 1 .36) 
1 .09 (0.87 to 1 .37) 


0.92 (0.69 to 1 .22) 
1.06 (0.83 to 1.36) 
1.17 (0.99 to 1.39) 
1.52 (0.91 to 2.55) 
1.17 (0.80 to 1.70) 
1.19 (0.97 to 1.47) 


0.73* (0.56 to 0.95) 
0.83 (0.65 to 1 .05) 
1.11 (0.97 to 1.28) 
0.95 (0.68 to 1.32) 
1 .05 (0.76 to 1 .45) 
1.13 (0.94 to 1.37) 


0.62** (0.44 to 0.87) 
0.71* (0.52 to 0.97) 
1.59** (1.18 to 2.14) 
0.81 (0.56 to 1.17) 
0.94 (0.62 to 1 .45) 
1.11 (0.89 to 1.38) 


0.86 (0.65 to 1.14) 
0.96 (0.75 to 1 .24) 
1.02 (0.86 to 1.22) 
1.22 (0.71 to 2.09) 
1.11 (0.76 to 1.61) 
1.11 (0.90 to 1.37) 



Values are ORs (95% CI). 

Models are adjusted for gestational age, gestational age squared, birthweight z-score, use of antenatal steroids, gender, infant year of birth and deprivation. 
*p<0.05, **p<0.01 , ***p<0.001 . 

tHigh volume was defined as being in the top quartile of units by number of care days provided to infants born at <32 +6 weeks gestation. 
BPD, bronchopulmonary dysplasia; NEC, necrotising enterocolitis; PMA, postmenstrual age; ROP, retinopathy of prematurity. 



rable 3 Descriptive statistics for the sample of preterm infants born <32 +6 weeks gestation by designation of the neonatal unit nearest to maternal place of residence 

Unit level designation Unit volume* 







Nearest unit 




p Valued, 




Nearest unit 




p Valuet, 




Nearest unit 


non-tertiary 




controlling for 


Nearest unit 


non-high 




controlling 




tertiary level 


level 


p Valuet 


deprivation 


high volume 


volume 


p Valuef for deprivation 


n (%) 


7167 (34.9) 


13 387 (65.1) 






7357 (35.8) 


13 197 (64.2) 






Gestation (weeks), mean (SD) 


29.6 (2.4) 


29.7 (2.3) 


0.040 


0.418 


29.6 (2.4) 


29.6 (2.3) 


0.181 


0.526 


Birth weight (g), mean (SD) 


1377.4 (429.2) 


1394.2 (424.5) 


0.007 


0.262 


1376.7 (426.7) 


1394.8 (425.7) 


0.004 


0.111 


Received full or partial course of antenatal 


4703 (65.6) 


8953 (66.9) 


0.069 


0.584 


4749 (64.6) 


8907 (67.5) 


<0.001 


0.052 


steroids 


















Deprivation score bottom 10% 


1751 (24.4) 


1611 (12.0) 


<0.001 


NA 


1476 (20.1) 


1886 (14.3) 


<0.001 


NA 


Male 


3820 (53.3) 


7165 (53.5) 


0.761 


0.854 


3958 (53.8) 


7027 (53.3) 


0.447 


0.378 


Birth in hospital with tertiary-level unit 


4753 (88.4) 


2290 (22.5) 


<0.001 


<0.001 


3839 (69.5) 


3204 (31 .9) 


<0.001 


<0.001 


Birth in hospital with high-volume unit 


3703 (68.9) 


3374 (33.1) 


<0.001 


<0.001 


4764 (86.3) 


2313 (23.0) 


<0.001 


<0.001 



All values are n (%) and are a proportion of the column total unless otherwise stated. 

*High volume was defined as being in the top quartile of units by number of care days provided to infants born at <32 +6 weeks gestation. 
tContinuous variables were tested by t test, categorical variables by x 2 test. 

i:p Value of F-test of coefficient on instrument from a regression of variable of interest on instrument and deprivation indicator. 
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infants with any missing data were excluded from the 
analyses (see online supplementary table Bl). 

The results remained robust to the inclusion of unit 
level random effects. We further excluded infants who 
died from analyses of the morbidity outcomes. This did 
not reveal any evidence of differences in the ORs except 
for the OR of treatment for ROP for infants admitted to 
tertiary-level care at the hospital of birth (OR 1.96, 95% 
CI 1.15 to 3.32, p=0.013; see online supplementary table 
B2). No evidence of an effect for the outcome defined 
as any in-hospital mortality and/or BPD was observed 
(see online supplementary table B2). Three alternative 
measures of volume were also used. In these sensitivity 
analyses, the OR of any in-hospital mortality remained 
significantly lower for very preterm infants admitted to a 
high-volume unit at the hospital of birth (see online sup- 
plementary tables B3 and B4). Only eight hospitals 
(4.8%) met the criteria of at least 100 VLBW infants per 
annum in any of the study years so that only a small pro- 
portion (6.5%) of the sample was inborn and admitted 
to these units. There is therefore imprecision around 
these results with wide CIs; among these infants, the OR 
of any in-hospital mortality was significantly lower but 
not statistically significant (see online supplementary 
table B4). 



DISCUSSION 

We examined the effects of designation and volume of 
neonatal care provided at the hospital of birth on mor- 
tality and morbidity outcomes for very preterm infants 
in England. Our key finding was a consistent reduction 
in the OR of mortality for very preterm infants admitted 
to high-volume neonatal units. We examined infants 
born at <26 +6 weeks gestation and those born at 27 +6 - 
32 +6 weeks gestation separately to reflect transfer policies 
and found a statistically significant reduction in the OR 
of mortality in the former group only. Furthermore, we 
found differences in the OR of mortality outcomes 
between standard logistic regressions and our preferred 
instrumental variables approach. The standard logistic 
regressions were generally found to underestimate the 
beneficial effects of high-volume care on mortality out- 
comes. This was expected given the aim of MCNs to 
transfer high-risk infants to high-volume and designation 
units. With regards to morbidity outcomes, treatment for 
ROP was the only morbidity for which a statistically sig- 
nificant effect was observed across analyses. We found 
that infants born at 27 +6 -32 +6 weeks gestation in hospi- 
tals with tertiary-level units were at increased OR of 
receiving treatment for ROP; however, only a very small 
number of these infants received treatment for ROP 
(86/17 995; 0.5%), suggesting the observed difference 
may not be clinically significant. 

Our preferred instrumental variables methodology, in 
the absence of a randomised assignment of infants to 
units, enabled us to estimate the causal effects of desig- 
nation and volume of neonatal care provided at the 
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hospital of birth using observational data. This approach 
has been widely applied in other healthcare evalua- 
tions. 26 However, we can only identify one previous 
application of this methodology to the evaluation of 
perinatal outcomes. 7 Our findings agree with the find- 
ings of a US-based study that examined the separate 
effects of level and volume of neonatal care. 4 We also 
found a reduction in the OR of mortality when analysing 
the annual number of VLBW admissions of inborn 
infants — a measure frequently used in US studies of this 
nature. 2 

We acknowledge limitations to our study. First, instru- 
mental variables methodology only identifies the effect 
of an intervention or treatment for those individuals 
whose assignment to treatment is altered by the instru- 
mental variable. 27 We do not know the effects for infants 
who would always be born in hospitals with a high-level 
or high-volume neonatal unit despite the location of the 
mother's residence (see online supplementary appendix 
A). Nevertheless, we demonstrated the validity of our 
instruments in meeting the required assumptions. 
Importantly, the assumptions required for an instrumen- 
tal variables methodology are weaker than those 
required to support a 'standard' analysis, which requires 
that infants are randomly assigned to hospitals of birth; 
otherwise the estimated ORs will be biased. Second, due 
to data limitations we cannot control for the effects of 
care and risk of death in the delivery suite at the hospital 
of birth. However, high-volume delivery units have been 
shown to be associated with a reduced risk of neonatal 
mortality. 28 29 Since high-volume delivery units are often 
found in hospitals with high-volume neonatal care this 
would lead us to suspect that our analyses underestimate 
the benefits of birth in hospitals with high-volume neo- 
natal care. 

Third, we are unable to disentangle the effects of the 
unit at the place of birth and subsequent transfers on 
final outcomes. We therefore cannot assess whether 
increasing the provision of transfers attenuates the 
increased OR of mortality associated with birth in hospi- 
tals without high-volume neonatal care. While identifica- 
tion of acute neonatal transfers was possible from our 
data, identifying the effects of transfer on outcomes pre- 
sents a number of difficult statistical issues. However, we 
expect that, if transfers to high-volume units reduce the 
OR of mortality, our effects presented in this paper 
underestimate the benefit of birth in a hospital with 
high-level or high-volume neonatal care (see online sup- 
plementary appendix A for an extended discussion), 
although neonatal transport itself may have negative 
effects on infant health outcomes. 30 31 A final limitation 
is that a small number of neonatal units in England 
(n=8) across MCNs do not contribute data to the NNRD 
and/or participate in NESCOR The effect of also includ- 
ing data from these units on outcomes remains a topic 
for future enquiry. 

An intervention that increases the proportion of very 
preterm infants born in hospitals with high-volume 



neonatal units may involve increasing the proportion of 
in utero transfers. Transfers of women prior to delivery 
are generally preferable because they are believed to be 
safer and less expensive than postnatal transfers of vul- 
nerable infants. 32 However, a study in 2009 showed that 
almost one-half of all in utero transfer requests to the 
London Ambulance Service were unsuccessful for non- 
clinical reasons. 33 Furthermore, studies from other coun- 
tries, including Portugal, Finland and the USA, have 
shown that in more regionalised systems as many as 90- 
95% of very preterm or VLBW infants are born in hospi- 
tals with tertiary designation neonatal units. 10 34 35 The 
effects of transfers within different organisational struc- 
tures for neonatal care remains an important area for 
future research especially as the new English 
Operational Delivery Networks will supersede the peri- 
natal MCNs as part of the changes following the Health 
and Social Care Act (2012). 36 

In conclusion, instrumental variables methodology did 
not reveal evidence of a difference in mortality out- 
comes between very preterm infants admitted to either 
tertiary or non-tertiary neonatal care at the hospital of 
birth. However, we do provide evidence of reduced OR 
of mortality for very preterm infants admitted to high- 
volume neonatal units at delivery hospitals. The effect of 
volume on neonatal outcomes is an important consider- 
ation for policy makers deciding the optimal organisa- 
tion of neonatal specialist services. 
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